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Introduction 



My goal is to describe new states of matter, which were initially encountered in purely theoretical 
investigations within the mathematical formalism used in physics, and subsequently were looked 
for and sometimes discovered in actual experiments. But first I want to remind you of an earlier 
instance of such a development advancing physics: Dirac's discovery of anti-matter. 

Dirac was seeking a relativistic and quantum mechanical equation to describe spin 1/2 electrons- 
free electrons to begin with. He proposed a matrix equation, first order in derivatives. Here it is. 



The Dirac equation is a beautiful equation of mathematical physics. It is also a beautiful 
equation in mathematics because, as Dirac himself stated, it provides a square root for the Laplacian 
of the Schrodinger equation, and the Laplacian is a much studied mathematical operator. 

But square roots come with both signs, and the energy eigenvalue E of the Dirac operator 
comes out positive, which is good, but also negative, which requires reinterpretation, because 
physical electrons carry only positive energy. 

With characteristic faith in the mathematical formalism, Dirac did not reject his equation 
with its negative energy modes. Rather he interpreted them as describing oppositely charged 
anti-electrons. This speculation was soon confirmed by the experimental discovery of positrons. 
Evidently keeping faith with a beautiful equation, more specifically, taking seriously the eigenvalue 
spectrum, no matter what its peculiarities may be, leads to physical discovery! 

There is no more news about Dirac's original equation. Current research proceeds because 
his equation has been extended and deformed in various ways. With these extensions we have 
encountered further spectral peculiarities, which are mathematically fascinating and whose physical 
interpretation exposes novel excitations in physical systems. 

Deformations of Dirac's equation arise from two motivations: mathematical and physical. Math- 
ematicians facing a beautiful equation instinctively look for the more general case, for an extended 
structure. On the other hand, physicists have established that Dirac-like equations describe low- 
energy excitations in actual condensed matter systems. Such descriptions are of course approximate, 



Dirac Matrix Equation 




(a • p + (3m) ^ = i f t ^ 

complex (charged excitations) 

a, (3 => matrices, p =^ i V, m =^mass 

energy eigenvalues: 

(a • p + fim) = E^ E 

E > empty in vacuum 



E < filled in vacuum 



m produces gap 

vacuum charge vanishes 
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but this does not detract from their usefulness in a restricted context. After all, even the original 
Dirac equation provides only an approximate treatment of relativistic electrons. 

Electron motion within materials extends over the volume occupied by the sample. However in 
certain materials electron motion is effectively confined to a line or to a plane. Correspondingly our 
Dirac equation needs to be reformulated on a one-dimensional line or on a two-dimensional plane, 
in addition to the three-dimensional volume. 

Even before delving into details, we can appreciate how aptly the Dirac equation confronts 
condensed matter phenomena: its positive and negative energy eigenvalues model free conduction 
electrons and bound valence electrons respectively. In the ground state the positive energy modes 
are empty because there are no conduction electrons, while the negative energy modes are filled with 
valence electrons. The mass term in the equation produces by its absolute value a gap between 
the conduction and valence bands, as in an insulator; a vanishing mass term leads to a gapless 
spectrum, as in a conductor. 

A profound deformation of the Dirac equation is achieved when the constant in space mass term 
is replaced by an inhomogenous term that varies with position. Only profiles that are significantly 
different from the homogenous case are interesting because small variations give rise to small, 
insignificant changes. To guide us in selecting interesting profiles, we classify them into distinct 
classes. Such classification uses concepts of topology and geometry and sparks the interest of 
mathematicians, while physicists are challenged to give a physical interpretation of the solutions 
to the Dirac equation in the presence of topologically non-trivial, position dependent mass terms. 

Examples of interesting profiles that can replace the homogenous mass term are: kink in one 
dimension, vortex in two and magnetic monopole in three. 



Dirac Equations 



[a ■ p + /3m]* E = E^e 

Dirac equation (mass term position-independent, homogenous) 
Usual: continuum solutions E > and E < 



Dirac equation in the presence of a defect 



(mass term position-dependent) 
[a -p + /3 </>(r)]* E = Eip E 



m — > (j)(r): phonon field 



Examples of Defects 



3-d magnetic monopole 



2-d vortex 



1-d kink-soliton 




Continuum solutions E > and E < and isolated E = solution 
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With Rebbi and Rossi we studied Dirac equations in these topologically nontrivial backgrounds. 
We found that in all instances the spectrum contains the expected positive and negative energy 
modes, but additionally there are also isolated zero-energy modes. Moreover the zero-energy so- 
lutions do not require a specific form for the inhomogenous mass profile; all that matters is that 
it belongs to a non-trivial topological class. It soon became clear that the presence of the isolated 
zero-energy modes can be established a priori by mathematical index theorems, which relate the 
occurrence of these modes to the topology of the mass term profile and the geometry of the space 
on which the equation is stated. 

This realization delighted us, and encouraged our belief that the mathematical elegance of the 
zero-mode solution ensures its profound physical significance. The remaining task was to describe 
the physical role of the zero modes. 



Polyacetylene 

The first and oldest instance of this story concerns polyacetylene. This is a lineal array of carbon 
atoms, separated by about one angstrom. The dispersion law connecting energy E to momentum 
p is sketched in the figure. The dashed line indicates the Fermi energy. 



1A | 1A 
• * • 



Fermi level 



Hopping in 1 dimensional array 

I E ( p ) linearization at Fermi level 
2 "Dirac Points" 




=^2x2 Dirac Hamiltonian kinetic term 
in 1 spatial dimension 



The energy curve intersects the Fermi level in two points, called Dirac points. Near these points 
the dispersion can be linearized. Thus an effective Hamiltonian is a 2 x 2 matrix multiplied by p, 
this provides the kinetic term of a one-dimensional Dirac Hamiltonian, but is only half the story. 

It is known that in polyacetylene the uniform array of carbon atoms is unstable. In fact 
equilibrium is achieved only after the atoms shift by about .04 angstroms to the left or to the right. 
Both cases are possible, giving rise to two equivalent domain structures. This is called the Peierls 
instability, and is modeled in our effective Dirac Hamiltonian by a homogenous mass term that 
takes either a positive or negative value. 



1A | 1A 

Unstable • • • f 



.04A 

-M 

Two stable b • • « « 



positions : 04A 



equilibrium 



With these steps we arrive at the Dirac-like equation for polyacetylene. 



4 



Polyacetylene Realization of 1-d Dirac equation 



cr 2 pMs E + a X (p^ E = EW E (2x2) 



linearization at 
Fermi level 



Peierls' instability 

constant m => mass gap \m\ 
b soliton 4> s zero mode U/ 



Pauli matrices 



When the mass term is inhomogenous we call it a phonon field, whose non-vanishing value 
describes the carbon shifts due to the Peierls instability. 

The energetics and the profiles of the phonon field are pictured in the figure. 

Energetics of Polyacetylene Phonon Field 




Energy density V((f>), as a function of a constant phonon field 0. The symmetric stationary point, = 0, is 
unstable. Stable configurations are at = + |m|, (A) and0 = — |to|, (B). 



Profiles of Polyacetylene Phonon Field 

(X) 



A 



B 




The two constant fields, ± | to |, correspond to the two domains (A and B). The kink-soliton fields, ±0 S , 
interpolate between the vacua and represent domain walls. 

The role of a kink-shaped mass profile is now apparent: it describes a domain wall, interpo- 
lating between positive and negative mass terms, that is between the two equivalent domains of 
polyacetylene. 
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A simple calculation shows that no matter what the detailed kink shape might be, the Dirac 
equation in the presence of a kink always supports a normalizable, zero-energy state. But in fact 
no explicit calculation is needed because an index theorem guarantees the result. 

The question remains how to treat the zero-energy state in the construction of the ground state 
of the material. Should the state be empty, like the positive energy states, or does it belong to the 
filled, negative energy states? Dirac has nothing to say about this, because he did not know about 
zero-energy modes! 

Fractional Charge Q 

Rcbbi and I gave the answer. Because there is no energy difference between the two possibilities, 
both have to be taken into account: empty and filled. Moreover the charge Q (or fermion number) 
of the empty state is —1/2, while the filled state carries charge +1/2. The fractional number is not 
an expectation value, but a sharp eigenvalue, that is there are no fluctuations. 

There are many ways to prove this surprising result. An intuitive argument is based on the fact 
that our effective Hamiltonian possesses a conjugation symmetry, such that the energy spectrum 
of filled states is the same as that of empty states, but the charge is opposite. Because of this 
symmetry our zero energy state lies in the middle of the gap, and the charge difference between 
empty and filled states is one. Since the charges should also differ only in sign they must take the 
values ±1/2. 

In a more formal argument, we construct the charge operator and find 
Q = contribution of E / states + \ (a^a — aa^) 
Here a) and a are the creation and annihilation operators for the zero mode, obeying the usual 
fermionic anti-commutation relations, acting on states to fill or empty them. 

a) a + aat =1 aa = 

| empty) = | filled) 
a | filled) = | empty) 
at (filled) = 
a | empty) = 
Q | state) = ±i | state) 

It follows that Q = ±1/2. 

Experimental confirmation came when Schrieffer and collaborators observed that charge frac- 
tionalization explains certain novel features of polyacetylene. Thereby they established the physical 
realty of this unexpected phenomenon. 

As a result physicists opened their minds to the occurrence of charge fractions. Thus when the 
quantum Hall effect emerged as a fascinating planar quantum system, fractionally charged states 
helped understand the relevant physics. 

But these fractions arise from a very different mechanism. There is no Dirac equation; there 
are no topological defects; most significantly, the quantum Hall effect exists in the presence of a 
magnetic field, which breaks time-inversion invariance , while our Dirac equations preserve time- 
inversion. 

Therefore, the question arises: beyond one dimension are there physical systems described by 
an effective Dirac equation that leads to charge fractionalization in a time-inversion symmetric 
manner. 
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Especially challenging is the fact that in higher dimensions the Fermi level defines a surface and 
its intersection with the energy dispersion is an extended manifold, unlikely to reduce to a finite 
number of Dirac intersection points. 

But then comes graphene. 

Graphene 

Graphene is a planar array of carbon atoms in a hexagonal lattice. It is a one-atom thick sample 
of graphite. The hexagon can be displayed as a superposition of two triangular sublattices, each 
possessing two Dirac points. 

Graphene hexagonal lattice 




Linearization at Fermi level 
2 "Dirac points" per sublattice 

Wallace showed this in 1947. He was a Canadian physicist, Canada's contribution to the 
American atomic bomb project, where he studied properties of graphite — an important neutron 
moderator. After the war work Wallace published his investigation of monolayer graphite, where 
he demonstrated that due to the specific band structure the energy momentum dispersion can be 
linearized near two Dirac points, just as in polyacetylene. Almost forty years later, Semenoff - 
another Canadian theorist — showed that the linear approximation leads to a Dirac kinetic term. 

The Dirac structure gives rise to specific experimental effects. These have been observed in 
graphene, and the 2010 Nobel Prize was awarded to the experimentalists Geim and Novoselov for 
their graphene work. 

But this is only half the way to fractional charge. Still needed is a mass-like term, which can 
take a constant value or an inhomogenous profile in the shape of a vortex. The constant value 
would produce a gap in the energy spectrum, while the topologically non-trivial configuration of a 
vortex would support a zero-energy mode. In other words, what is needed is a planar analog of the 
Peierls instability. 

Chamon and collaborators have identified for graphene the possibility of such an effect, called a 
Kekule structure. It preserves time-inversion invariance and bonds carbon atoms in the plane with 
a pattern that generalizes the one-dimensional Peierls instability. 
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Graphene hexagonal lattice with Kekule distortion 





But unlike the Peierls effect there is no guarantee that the Kekule structure is always present. 
If it occurs then the effective Dirac equation reads 

Graphene Realization of (2-d) Dirac Equation 



cr ■ p 

A* 



A-^ 

<7 • p 



E 



(4x4) 



linearization 
at Fermi level 



Kekule distortion 
A constant mass gap 
A soliton/vortex =>■ mid-gap state 



With a vortex profile for A , Chamon used the solution to a related equation, which was found 
earlier by Rossi and me, to display a zero-energy mode that leads to fractional charge, in the 
same manner as for one-dimensional polyacetylene. Thus in spite of the tenuousness of the Kekule 
structure, we have here at the very least an existence proof that fractional charge can also arise in 
the plane in a time-inversion symmetric scenario. 

Major ana Fermions 

Complex fields describe charged particles and their anti-particles of opposite charge. This is as in 
the Dirac equation for electrons and positrons. In nature there are also neutral particles that are 
their own anti-particles. Examples are 

neutral pion (spin 0) 
photon (spin 1) 
graviton (spin 2) 

All these are integer-spin bosons, and real fields describe them. 
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The question arises: can there exist neutral half- integer spin fermions that are their own anti- 
particles? Thus far experimentalists have not identified any, but theorists like them very much: 
they use them in speculations on particle physics where they are supersymmetric partners of neutral 
bosons; in cosmology they possibly constitute dark matter and describe neutrinos; most recently 
they have arisen in condensed matter physics as exotic superconducting states. 

Neutral, spin 1/2 particles that are their own anti-particles are called "Majorana fermions" 
because the mysterious Italian physicist Majorana discovered the equation that describes them. 
Let me review. 

Majorana Equation 





Majorana Matrix Equation 
(a • p + (3m)V = i§t X S> 
^ real (neutral excitations) 
p = i V imaginary 

% Imaginary } Majorana Representation 



Because they are spin 1/2 particles, the Dirac equation remains valid. Majorana realized that 
there exists a representation of the Dirac matrices a and f3 , where a. is real and (3 is imaginary. 
In this representation, called the "Majorana representation," one can demand that the spinor ^ 
is real, and then the Dirac equation becomes the real Majorana equation. This ensures charge 
neutrality of a spin 1/2 particle, identified with its anti-particle. 

The Majorana equation is a further step in the deconstruction of a differential operator: real 
second order Laplacian; complex first order Dirac operator; real first order Majorana operator. 
Thus it possesses mathematical interest in addition to its physical role 

In fact matrices arising in the Majorana equation need not possess the above-stated reality 
properties. Since similarity transformations on matrices do not affect physical content, it is sufficient 
to demand in any representation that there exists a matrix C such that 



Ccx*C- 1 = cx, CI3*C~ 



-P 



and the reality condition on ^ is replaced by the pseudo-reality requirement C^S>* = ^fr. In the 

-- 2 



Majorana representation C = I, while with the usual, complex a and /3, C=\ . 2 



■ur 

A four-component spinor satisfying the pseudo-reality condition necessarily takes the form 



and the four-component equation reduces to the two component Majorana equa- 
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tion. 

(j ■ p + i a 2 ip ib* = i — ib 
at 

Note that the mass term couples ib to ib*. Therefore the equation does not entail charge or 
particle conservation; indeed excitations carry no charge. 

Recent developments in neutrino particle physics have established that neutrinos have mass. 
The occurrence of neutrino oscillations indicates that individual species of neutrinos are not con- 
served. This suggests that neutrinos are Majorana fermions. Nowdays there is intense search for 
definite verification of the Majorana nature of the neutrino, making the Majorana equation relevant 
to physics in three dimensions. 

The Majorana equation taken in two spatial dimensions allows replacing the constant mass 
term with an inhomogenous vortex profile, whereupon the equation possesses zero modes, as was 
first found by Rossi and me. 

Remarkably this same planar equation has recently been posited by Fu and Kane in their 
description of a superconductor in contact with a topological insulator. A proximity effect causes 
Cooper pairs to tunnel through the surface of the topological insulator. This effect is modeled by 
the planar Majorana equation, which also supports the above mentioned zero mode in the presence 
of the vortex. 

Owing to its Majorana nature, the zero mode does not carry any charge. Moreover, it entails 
a novel mode algebra for its mode operators, like the ones introduced earlier in the polyacetylene 
story. 

aa^ + a^a = 1 but a = a^ (Majorana) 
=>- aa = \ 

This novel algebra is neither conventionally fermionic nor bosonic. Its action calls for two states, 

l+>,i->. 

«i+> = ^i-> °i-> = ^i+> 

It has been suggested that this algebra will find application in quantum computation. This is 
because the zero mode is stable against perturbations, since topological considerations guarantee 
its presence. 

Now we wait to learn who will first discover Majorana fermions; particle physicists with neutrino- 
less double beta decay for example, or cosmologists by identifying dark matter, or condensed matter 
people with exotic superconductors, or will quantum computers beat out the physicists. 

Into the Future 

What more is there to do? We have discussed one-dimensional kinks in polyacetylene, two- 
dimensional vortices in graphene and in exotic superconductors. What about three dimensions, 
where Dirac originally posited his equation? 

With Rebbi, we showed that replacing the constant mass term by a three-dimensional profile 
of a magnetic monopole again leads to zero-modes in the extended Dirac equation. Since charge 
remains conserved, this system exhibits fractional charge. 

Unfortunately we do not have a good physical candidate for a magnetic monopole, even though 
material scientists have occasionally discussed effective monopoles, like the ones in spin-ice. But 
these are not topological defects. 
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Therefore a good future problem would be to determine whether magnetic monopoles exist and 
whether they induce fractional charge in three dimensions. 

This work is supported by the U.S. Department of Energy under cooperative research agreement 
Contract Number DE-FG02-05ER41360. 



11 



